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Abstract 

This paper is the second in a series revisiting the (effect of) Faraday rotation. We formulate 
and prove the thermodynamic limit for the transverse electric conductivity of Bloch electrons, 
as well as for the Verdet constant. The main mathematical tool is a regularized magnetic and 
geometric perturbation theory combined with elliptic regularity and Agmon-Combes- Thomas 
uniform exponential decay estimates. 
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1 Introduction 



1.1 Generalities 

The rotation of the polarization of a plane-polarized electromagnetic wave passing through a mate- 
rial immersed in a homogeneous magnetic field oriented parallel to the direction of propagation, is 
known in physics as the Faraday effect (called sometimes also Faraday rotation) . The experiment 
consists in sending a monochromatic light wave, parallel to the Oz direction and linearly polarized 
in the plane xOz. When the light enters the sample, the polarization plane starts rotating. A 
simple argument based on (classical) Maxwell equations shows that there exists a linear relation 
between the angle 9 of rotation of the plane of polarization per unit length and the transverse 
component of the conductivity tensor of the material (see e.g. formula (1) in 29J). For most ma- 
terials - and we will restrict ourselves to this case - the transverse component of the conductivity 
tensor vanishes when the magnetic field is absent and is no longer zero when the magnetic field is 
turned on. Under the proviso that the dependence of the conductivity tensor upon the strength 
B of the magnetic field is smooth, for weak fields one expands the conductivity tensor to the first 
order and neglect the higher terms. The coefficient of the linear term is known as the Verdet 
constant of the corresponding material. 

It follows that the basic object is the conductivity tensor and the main goal of the theory 
(classical or quantum) is to provide a workable formula for it, in particular for the Verdet constant. 
The problem has a long and distinguished history in solid state physics theory and the spectrum 
of possible applications ranges from astrophysics to optics and general quantum mechanics (see 
e.g. [29, 28, 33, 13, 23, and references therein). 

Using quantum theory in the setting in which the sample is modeled by a system of independent 
electrons subjected to a periodic electric potential, Laura Roth [29] obtained (albeit only at a 
formal level) a formula for the Verdet constant in full generality and applied it to metals as well as 
semiconductors. Roth's method is based on an effective Hamiltonian approach for Bloch electrons 
in the presence of a weak constant magnetic field (see |31] and references therein) which in turn 
is based on a (proto) magnetic pseudodifferential calculus (for recent mathematical developments 
see [21] and references therein). 

But Roth's theory is far from being free of difficulties. Due to her highly formal way of doing 
computations, it seems almost hopeless - even with present day mathematical tools - to control 
the errors or push the computations to higher orders in B except maybe the case of simple bands. 
Even more, the final formula contains terms which are singular at the crossings of Bloch bands. 
Accordingly, in spite of the fact that it has been considered a landmark of the subject, it came as 
no surprise that at the practical level this theory only met a moderate success and a multitude of 
unrelated, simplified models have been tailored for specific cases. 

Our paper is the second in a series aiming at a complete, unitary and mathematically sound 
theory of Faraday effect having the same generality as Roth's theory (i.e. a theory of the conduc- 
tivity tensor for electrons subjected to a periodic electric potential and to a constant magnetic 
field in the linear response approximation), but free of its shortcomings. More precisely in the first 
paper [13j we started by a rigorous derivation of the transverse component of the conductivity 
tensor in the linear response regime for a finite sample. It is given as the formula 11.101 below . To 
proceed further, we employ a method going back at least Sondheimer and Wilson [32] and which 
has been also used in the rigorous study of the Landau magnetism [TJ [Tni [5l [H [7] ■ The basic idea 
is that the traces involved in computing various physical quantities can be written as integrals 
involving Green functions (i.e. integral kernels in the configuration space of either the resolvent 
or the semi-group of the Hamiltonian of the system) , which are more robust and easier to control. 

As it stands, the conductivity tensor depends upon the shape of the sample and of boundary 
conditions which define the quantum Hamiltonian. The physical idea of the thermodynamic limit 
for an intensive physical quantity is that in the limit of large samples it approaches a limit which 
is independent upon the shape of the sample, boundary conditions etc. The existence of the 
thermodynamic limit is one of the basic (and far from trivial) problem of statistical mechanics 
(see e.g. [SO] [3]). In [13] we took for granted that the thermodynamic limit of the the transverse 
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component of the conductivity tensor exists and the limit is smooth as a function of the magnetic 
field strength B. Moreover, we assumed that the thermodynamic limit commuted with taking the 
derivative with respect to B. Under these assumptions, we gave - among other things - explicit 
formulas for the Verdet constant in terms of zero magnetic field Green functions, free of any 
divergences. 

The proof of the thermodynamic limit, which from the mathematical point of view is the most 
delicate part of the theory of the Faraday effect, was left aside in [13j, and is the content of the 
current paper. The mathematical problem behind it is hard due to the singularity induced by the 
long range magnetic perturbation. Even for a simpler problem involving constant magnetic field 
- namely the Landau diamagnetism of free electrons - the existence of the thermodynamic limit 
leading to a correct thermodynamic behavior was a long standing problem. Naive computations 
led to unphysical and contradictory results (see [T] for historical remarks). Accordingly, the first 
rigorous results came as late as 1975 [T] and were based on various identities expressing the gauge 
invariance which was crucial in dealing with the singular terms appearing in the thermodynamic 
limit. Even though the importance of gauge invariance was already highlighted in [Tj, an efficient 
way to implement this idea at a technical level was still lacking. Only recently a regularized mag- 
netic perturbation theory based on factorizing the (singular in the thermodynamic limit) magnetic 
phase factor has been fully developed in [10l[TTl[l2, 27 . This regularized magnetic perturbation 
theory has been already used in [101 O El E] in order to prove far reaching generalizations of the 
results in [Tj. 

Coming back to the Faraday effect, we would like to stress that the object at hand is much more 
singular than the one encountered in the Landau diamagnetism. This adds an order of magnitude 
to the mathematical difficulty and requires an elaborate and tedious combination of regularized 
magnetic perturbation theory with techniques like Combes-Thomas exponential decay, trace norm 
estimates and elliptic regularity. 

We expect that the method developed here in order to control the thermodynamic limit in the 
presence of an extended magnetic field to be useful in related problems, e.g. to obtain an elegant 
and complete study of the diamagnetism and de Haas-van Alphen effect for electrons in metals. 

The content of the paper is as follows. In the rest of this Introduction we state the mathematical 
problem, give the main result in Theorem 1 1.11 and since the proof is quite long and technical we 
will briefly describe the main points. The other sections are devoted to the proof of our main 
theorem. The core of the proof heavily involving regularized magnetic perturbation theory is 
contained in Section 4. Some technical estimates about exponential decay with uniform control in 
the spectral parameter are given as an Appendix. 

1.2 The main result 

Consider a simply connected open and bounded set Ai C K"^, which contains the origin. We 
assume that the boundary dAi is smooth. Consider a family of scaled domains 

AL^{xeR^ : x/LeAi}, L>1. (1.1) 

We have the estimates 

Vol(Ai) - L^, Area(aAL) - L^. (1.2) 

We are interested in the thermodynamic limit, which will mean L oo, that is A^ will fill out 
the whole space. The one particle Hilbert space is Hl '■— L^{Al)- Note that we include the case 
L — oo. 

The one body Hamiltonian of a non-confined particle, subjected to a constant magnetic field 
(0, 0, B), in an external potential V, formally looks like this: 

i/oo(s) = p2(B) + y, (1.3) 

with 

F{B) = -iV + Ba. = P{0) + BsL. (1.4) 
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Let us explain the various terms. Here a(x) is a smooth magnetic vector potential which generates 
a magnetic field of intensity B — 1 i.e. V A a(x) = (0,0, 1). A frequently used magnetic vector 
potential is the symmetric gauge: 

A(x) = ina Ax= (-a;2/2,a.i/2,0), (1.5) 

where is the unit vector along z axis. We neglect the spin structure since it only complicates 
the notation and does not influence the mathematical problem. 
On components, (|1.4p reads as: 

Pj{B) = Dj + Baj =: Pj{0) + Ba,, j G {1, 2, 3}. (1.6) 



We will from now on assume that is a C°°(K'^) function, periodic with respect to the lattice 
Z^. Standard arguments then show that Hao{B) is essentially self-adjoint on C^(R^). 

When L < oo we need to specify a boundary condition. We will only consider Dirichlet 
boundary conditions, that is we start with the same expression as in ()1.3|) . defined on C^{Kl), 
and we define Hl{B) to be the Friedrichs extension of it. This is indeed possible, because our 
operator can be written as — Ajj + W, where A^j is the Dirichlet Laplacian and is a first order 
differential operator, relatively bounded to —Ad (remember that L < oo). The form domain 
of Hl{B) is the Sobolev space HI{Kl), while the operator domain (use the estimates in section 
10.5, Lemma 10.5.1, in [19]) is H'^{Kl) n HI{Kl). Moreover, Hl{B) is essentially self-adjoint 
on C^^^(Ai), i.e. functions with support in A^ and indefinitely differentiable in A^ up to the 
boundary. 

Another important operator is (— iV + i?a)|,, i.e. the usual free magnetic Schrodinger operator 
defined with Dirichlet boundary conditions. We know that its spectrum is non-negative for all 
L > 1. By adding a positive constant, we can always assume that the spectrum of Hi^{B) is 
non- negative, uniformly in L > 1. 

Let us now introduce the physical quantity we want to study. Consider w e C and 3(w) < 0. 
For some fixed /i € M and /3 > 0, define the Fermi-Dirac function on its maximal domain of 
analyticity: 

Define 

. ( |Imw| 1 
rf:=mm|-,^|, (1.8) 

and introduce a counter-clockwise oriented contour given by 

V^^ = {x ±id : a < X < oo} {a -f : —d < y < d} (1-9) 

where a + 1 lies below the spectrum of Hl{B). By adding a positive constant to V , we can take 
a — —1 uniformly in L > 1 and B g [0, 1]. 

We introduce the transverse component of the conductivity tensor (see [IHIIIS]) as 

'^-(^^ = - Volk) ^'-''^ 
■ TtJ^ fpoiz) {Pi{B){Hl{B) - z)-^P2{B){Hl{B) - z - to)-' 

+ Pi{B){Hl{B) - z + uj)-^P2{B){Hl{B) - z)-^] dz. 

Here Tr assumes that the integral is a trace-class operator. Now we are prepared to formulate 
our main result. 

Theorem 1.1. The above defined transverse component of the conductivity tensor admits the 
thermodynamic limit; more precisely: 
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i. The following operator, defined by a B(L'^{Al))- norm convergent Riemann integral, 



Fl := ^ fFDiz){Pi{B){HLiB) - z)-^ P2{B){Hl{B) - z - uj)'^ 

+ Pi{B){Hl{B) - z + uj)-^P2{B){Hl{B) - z)-^}dz, (1.11) 
is in fact trace-class; 

ii. Consider the operator F^o defined by the same integral but with H^oiB) instead of H]^(B), 
and defined on the whole space. Then is an integral operator, with a kernel J^(x, x') which 
is jointly continuous on its variables. Moreover, the continuous function defined by M.'^ 3 x. ^ 
sb(x) := J^(x, x) e R is periodic with respect to 1? ; 

iii. Denote by $7 the unit cube in R^. The thermodynamic limit exists: 

a^[B) := lim ^^(B) = - / .SB(x)dx. (1.12) 
Moreover, the mapping B ^ sb ^ L°°{fl) is differentiable at B = and: 



9b Coo (0) = - / dsSB 



(x)dx = lim dB(TL{0)- (1.13) 



B=0 ^ 



Remark 1. The formula (|1.12p is only the starting point in the study of the Faraday rotation. A 
related problem is the diamagnetism of Bloch electrons, where the main object is the integrated 
density of states of magnetic Schrodinger operators (see [171 [181 [20] ) . For a systematic treatment 
of magnetic pseudo-differential operators which generalizes our magnetic perturbation theory, see 

[HI [M [211115]. 

Remark 2. The Dirichlet boundary conditions are important for us. Even though we suspect 
that our main result should also hold for Neumann conditions and for less regular domains (see 
[51[T3]), we do not see an easy way to prove it. 

Remark 3. We believe that the method we use in the proof of (|1.13p can be used in order to 
obtain a stronger result: the mapping B ^ sb & L°°{il) is smooth and for any n > 1: 



d^a^{B) = - f 9]^ss(x)dx = lim d^aUB). 



(1.14) 



We leave this statement as an open problem. In the rest of the paper we give the proof of Theorem 

[TH 



1.3 A short description of the proof strategy 

Since the proof is rather long, we list here the main steps and ideas. Let us start with some general 
considerations about the thermodynamic limit. 

If we are interested in the thermodynamic limit of a quantum physical quantity, the object 
we need to control is the trace of the operator representing the corresponding quantity. The 
basic ideea consists in writing this trace as an integral of the diagonal value of the operator's 
integral kernel (Schwartz kernel) over the confining box. This procedure makes the quantum 
thermodynamic limit look very similar to what happens in classical statistical mechanics. More 
precisely, we need to show that the difference between the integral kernel for the finite box and 
the one for the entire space decays sufficiently rapidly with the distance from the boundary of the 
box, so that the replacement of the integral kernel for the box with the one corresponding to the 
entire space gives an error term increasing slower than the volume, which then disappears in the 
limit. 

It turns out that for the transverse conductivity this kernel is far more complicated than say the 
heat kernel - whose behavior has been extensively studied in the literature. Notice for example, that 
the integrand in (jl.lOp contains two resolvents sandwiched with magnetic momentum operators. 
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Thus we need a good control of their integral kernels, in particular when the distance between 
their arguments increases to infinity, and all that uniformly in the spectral parameter z. Since a 
constant magnetic field is present, the biggest difficulty is to deal with the linear growth of the 
vector potential. Here, the use of gauge covariance is crucial. 

Now let us list the main ideas of the proof. For the first statement of the theorem one simply 
uses integration by parts with respect to z in order to transform the integrand into a product of 
Hilbert-Schmidt operators. 

The proof of the second statement is based on elliptic regularity. The main technical difficulty is 
to control the z behavior of all our bounds, especially the exponential localization of the magnetic 
resolvents sandwiched with momentum operators. We also have to control the linear growth of 
the magnetic potential. We turn the operator norm bounds which we obtain in the Appendix into 
pointwise bounds for certain integral kernels. It is a long road using magnetic perturbation theory, 
but nevertheless, we use nothing more than well-known Combes-Thomas exponential bounds, local 
gauge covariance, the Cauchy-Schwarz inequality, and integration by parts. 

The third statement of Theorem 11.11 contains the main result and is proved in Section 4. We 
start with a bit strange three-layered partition of unity defined in (|4.1[) - (|4.12[) . This idea goes back 
at least to [H [5] . The main effort consists in isolating the bulk of - where only operators defined 
in the whole space will act - from the region close to the boundary. Note that in the absence of 
the magnetic field, it would be enough to work with only two cut-off functions: one isolating the 
bulk from the boundary, and the other one supported in a tubular neighborhood of the boundary. 
When long-range magnetic fields are present, this is not enough. The tubular neigborhood needs 
to be chopped up in many small pieces, in order to apply local gauge transformations (see below 
why we need them). 

The central idea in proving (|1.12p is to show that the contribution to the total trace of the region 
close to the boundary grows slower than the volume. Technically, this is obtained by approximating 
the true resolvent {Hl{B) — z)~^ with an operator Ul{B,z) given in (|4.16p . Ul{B,z) contains 
the bulk term, plus a boundary contribution which consists from a sum of terms each locally 
approximating {Hl{B) — z)~^ and containing a specially tailored local gauge given in (|4.13p . 
These locally defined vector potentials are made globally bounded with the help of our third layer 
of cut-off functions g-y's. The switch to the local gauge is performed through the central identity 
p.l8[) . An explanation of why Ul{B, z) is a convenient approximation for the full resolvent can be 
found right after (|4.18p . and there is the place where we first fully use the exponential estimates 
of the Appendix. 

In Proposition l4. ll we prove that we can replace {Hl{B)—z)^^ with Ul{B, z) in the conductivity 
formula, without changing the value of the limit. In Proposition 14.11 we show that only the bulk 
term from Ul{B^z) will contribute. In Proposition 14.61 we show that removing the cut-offs only 
gives a surface contribution. 

The proof of (|1.13p is heavily based on magnetic perturbation theory. Although the technical 
estimates are considerably more involved than at the previous point, the main idea is the same: the 
boundary terms can be discarded. The full power of the magnetic phases is used in Lemma [4.101 
an heuristic explanation of how and why they manage to kill the linear growth of the magnetic 
potential is given right after (|4.73p . 

2 Proof of i. 

The integrand in (jl.lip is a bounded operator, with an norm bounded by a constant times 
|5R(z)P, uniformly in L, as we can see from (|1.6p and (|5.2p . Because fpD has an exponential decay 
in 1 3^(2:) I, the integral converges and defines a bounded operator. Let us note that the integrand is 
not a trace class operator under our conditions. But the total integral is a different matter. The 
point is that we can integrate by parts with respect to z by using anti-derivatives of fpD which 
are still decaying exponentially at infinity. By doing this at least four times, we obtain integrals 
of the form 
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^ f{z)P,{B){HL{B)-z)-"^P,{B){HL{B)-z-u;)-"dz (2.1) 

where m + n > 5, hence max{m,n} > 3. Assume that m > 3. Then we can write the above 
integral as 

^ fiz)P,iB){HLiB)-z)-"'+'{HLiB)-z)-^P2iB)iHLiB)-z-u;)-"dz. (2.2) 

The main point is that (Hl{B) — z)^^ is Hilbert-Schmidt since we can write 
{HUB)-z)-' 

= ii-tV + Bb.)1 + {[(-zV + B2.)l + 1]{Hl{B) - z)-'] , (2.3) 

and by using (|5.7p with 6 = 0, and (|5.13p . we obtain 

WiHLiB) ~ z)-^\\b, < const • /VoI(iU)(3fi(z)), (2.4) 

where the above constant does not depend on L and z. Thus {Hl{B) — z)^"^ is trace class, and 
the trace norm of the integrand in (|2.2p is bounded by 

const - I /(z) I • (3?(z))* • VoI(Al) 

where again the constant is uniform in L and z. This now is integrable on the contour, thus the 
integral defines a trace class operator. Moreover, its trace grows at most like the volume of A^, 
hence limsup^^g^ IctlI < oo. □ 

Remark. The same type of argument may be used to show that <tl{B) is smooth in B, by 
repeatedly using the formal identity 

dB{HL{B) - z)-i = -{Hl{B) - z)-\dBHL{B)}{HL{B) - zy^ (2.5) 

in the sense of bounded operators. Note the important fact that dsHLiB) will generate some 
linear growing terms coming from the magnetic vector potential a(x), therefore the trace norm 
of the new integrand will grow like L*. We therefore cannot conclude here that the derivatives 
|9^(Tl(-B)| will admit a finite limsup when L — > oo. 

3 Proof of ii. 

We are going to prove the regularity statement for the kernel without using the periodicity of V , 
only the fact that the potential is smooth and bounded on R^^ together with all its derivatives. 

The strategy consists in integrating by parts with respect to z many times, such that we obtain 
high powers of the resolvent {Hoo{B)-C,)-^ . Then we wiU prove that Pj{B){Hoo{B)-C)-^ Pk{B) 
has a smooth kernel which does not grow too fast with (|C|)- 

Let us now be more precise and start with some technical results. 

Proposition 3.1. Fix < 77 < 1 and choose z e C with dist{z, [0, 00)} = rj > {]. Let r = (|3?(z)|). 
Then the operator {Hoo{Q) — has an integral kernel G'i(x, x';z) which is smooth away from 
the diagonal x = x'. There exists 6 > and some M > 1 such that for any multi-index a G 
with \a\ < 1 we have the estimate 

sup |x-x'|l"l+ie<^l"-"'l|D2Gi(x,x';z)| = Ci (a, r;)(r)*^ < 00. (3.1) 
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Proof. The result without the exponential decay is essentially contained in [T5]. The symbol of 
ffoo(O), denoted by /io(x, belongs to Sl^iR^ x R^), 7Joo(0) S i?,o(M^) (see Example 3.1 in [Tg), 
and is uniformly elliptic. 

Fix A > large enough. We can apply Theorem 4.1 in |15| and construct a parametrix for 
Hoo{0) + A starting from the symbol go(x,C) := l/(/i(x,0 + A) e •^^^^(R^ x R^). The symbol 
giving the parametrix is an asymptotic sum of symbols, starting with qq, then the next one is in 
and so on. Each term gives a contribution to the integral kernel of the parametrix. The most 
singular contribution is (in the sense of oscillatory integrals) : 

e^(''--')-«go(x,e)rfe 



{2nf 

We only have to consider a few terms besides this one, since symbols in S^q for large N generate 
more and more regular kernels at the diagonal. By standard "integration by parts" tricks, and 
using the fact that we work in three dimensions, one can prove the estimate 

sup |x-x'|l"l+i|L>2Gi(x,x';-A)| =const(a,A) < cx). (3.2) 

In fact, outside the region |x — x'| > 1 we can integrate by parts several times with respect to ^ 
and prove that Gi(x, x'; — A) decays faster than any power of (|x — x'|). But the Combes- Thomas 
method will give a better, exponential localization. 

The important thing is that the estimates from the Combes-Thomas argument can be 
transferred into point-wise estimates for the kernel. Let us now prove this. 

Using (j5.13p at L — oo and B = 0, together with the triangle and Cauchy-Schwarz inequalities, 
we get that (— A-|-A)^^ with exponential weights maps into L°° . The key estimate is (0 < c < 1) 

r- r- , „-(l-c)\/A|x-x'| 
-cVA|x-xo| I 7^ I „'^|„cVA|x -xo| <- 

e |Aoo(x,xj|e < ■ 

Since we can write: 

(Hoo(0) + A)-i = (-A + A)-i(-A + A)(i/oc(0) + A)-\ (3.3) 

using (j5.12|) (at L = oo and S = 0), it follows that this resolvent with exponential weights is a 
bounded map from into More precisely, for any xq € R^, there exists < c < 1 small 

enough such that: 

sup ||e-=^<-^°>(i7oo(0) + A)-ie^^<-'"'>||B(L2,Loo) < const, A > Aq. (3.4) 



Now if we look at the map 



c-( 



^) 9 * ^ / Gi(xo,x;-A)e*^<^-^°>*(x)dx, 

JR3 



(it makes sense to fix xq since the resolvent maps smooth functions into smooth functions), we 
see that by using (|3.4p we can extend this map to a linear and bounded functional on L? . Riesz' 
representation theorem then gives: 

sup ||e=^<-''«>Gi(xo,-;-A)|U2 = sup ||e^^<-^">Gi(-, xq; -A)|U2 < const, (3.5) 

xoeM3 xoSM^ 

uniformly in A > Aq. Using this, together with the Cauchy-Schwarz and the triangle inequality, 
we get that the integral kernel G2(x, x'; —A) of (i/oo(0) + A)~^ obeys uniformly in A > Aq: 

supe^^l''"'''l|G2(x,x';-A)| (3.6) 

x,x' 

< sup / |e'=^<''-''">Gi(x,x";-A)e'=^<^""^'>Gi(x",x';-A)|dx" < const. 

X,x' Jr3 
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Now if |x — x'l > 1, write 

G2(x,x';-Ai)dAi, (3.7) 

which together with (|3.6p and the integrabihty of e~5v^ imply that 

sup e^^l''-^'l|G'i(x,x';-A)| < const. (3.8) 

|x-x'|>l 

We can also deal with derivatives with respect to x, by showing that the operator Dj {Hoc (0) + 
A)"^ {N large enough) has an integral kernel DjGAr(x, x'; —A) obeying the same type of estimate 
as in (|3.6p . This is done by commuting Dj several times with a few resolvents; let us see how it 
works. 

First, by commuting we have: 

Dj{H^{Q) + A)-i = (i/oo(0) + Xy^Dj + (ff^(O) + A)-2Ti, (3.9) 
Ti = [H^{Q),Dj] - [ffoc^(O), [Hoo{0),D,]]{H^{0) + \)-\ 

where Ti is bounded due to the fact that [i7oo(0), -Dj] is bounded, while [iJoo(O), [i?oo(0), i'j]] is 
relatively bounded with respect to i?oo(0). 
Second, by commuting twice we have: 

DjiH^iO) + A)-2 = (i7^(o) + xy^Dj + (i/oo(0) + Xy'^Ti + (i?oo(0) + Xy^Ti{Hoo{0) + X)-\ 

(3.10) 

This identity allows us to write: 

Dj{H^{0) + X)-^+' = (iJoo(O) + A)-2r(iJ^(0) + A)-^+5, (3.11) 
T = {D, + (iJoo(O) + A)-iri + Ti(i/^(0) + A)-i}(i?oo(0) + A)-\ 

where T with exponential weights is bounded from to L^. Then we prove that the integral 
kernel of Dj{Hoo{0) + A)^^+^ obeys an estimate like in p.Sp . then from the identity 

oyH^O) + A)-^ = D,{H^{0) + xy^+^H^o) + A)-2 

we get the needed L°° estimate by mimicking (|3.6p . 
Then we write 

D,Gii^,x';-X)= \ \ / dXj dX2... dAjvi?jGAr(x,x';-A') 

and propagate the exponential decay over the integrals in A. 

Therefore we can state the first result regarding the exponential localization. For A large 
enough we have: 

sup |x-x'|l"l+iel''-'''l|i:i2Gi(x,x';-A)| =const(a,A) < oo. (3.12) 

Now let us investigate the z dependence. Let us apply the resolvent identity several times and 
get (iV > 2): 

{H^{Q) - zy^ = (H^iO) + A)-i + (z + A)(i/eo(0) + A)-^ + . . . 

+ (z + A)^(i7oo(0) + A)-^(i/oo(0) - zy\ (3.13) 

The idea is to keep the z dependence to the right in the last term, and to keep a regular kernel to 
the left. We start with a norm estimate. From the usual resolvent identity: 

(i7oo(0) - z)-i = (i/oo(0) + A)-i + (z + A)(i7oo(0) + xy\Hoo(.o) - zy\ 
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the use p.Sp and (|5.5p (with L = oo) provides us with some Ni > N such that: 

sup||e-^<-''«>(^oo(0) - zr'e^<-^°^\\BiL-,L^) < const(r;) • (3.14) 

xo 

This estimate impHes that the map (initially defined on compactly supported functions) 

L2(m3)3,j,^ f Gi(xo,x;z)e^<^-^o)^'(x)dx e C, 
Jr3 

is a bounded linear functional. Riesz' representation theorem leads us to: 

sup ||e^<-"°>Gi(xo,-;^)|U2 = sup ||e^<-"°>Gi(-, xq; z)|U2 < const(?7) • r^^ (3.15) 
xoeK^* xoeR3 

We are only left with the case in which we have a derivative on the left. Using (|3.15p in (|3.13p 
and the other results we have obtained for the kernel where z = A, it is not hard to obtain the 
exponential decay claimed in p.ip . □ 



Proposition 3.2. Assume that a £ {0,1} and < 77 < 1. Let dist{z, [0, 00)} — rj and N is 
large enough. Then the operator P^{B){H^{B) ~ z)-^ P^-''{B) has a jointly continuous integral 
kernel ii'Ar,B(x, x'; z), and there exists (5 > small enough and M large enough such that: 

sup e<^l''~'''l|ii:Ar,B(x,x';z)| < consi(A^,B,?7) • (r)*^ (3.16) 

x,x'GR3 

Proof. Although this particular result might be obtained with other methods, we will employ 
the magnetic perturbation theory as developed in [111 [5l [27] . For different approaches involving 
magnetic pseudo-differential calculus, see pT l ITS l [24 l [25 ] [21] . 

We can assume that the magnetic vector potential is expressed in the transverse gauge, given 
in (|1.5p . Define the antisymmetric magnetic gauge phase 

'^o(x,y) :== -A(y) • x = ]^{v2Xi - yixa) ^63 • (x Ay), (3.17) 

where 63 denotes the unit vector (0,0,1) £ M'^. Then we have the following identity (true for 
instance on Schwartz functions), valid for every vector y kept fixed in R^^: 

{Px(0) + SA(x)}e*'^'^"(^'y) = e^-^'^«('''yHPx(0) + BA(x - y)}. (3.18) 

For every z G C \ [0, 00) define 

5s(x, x'; z) := e'^'^°("'"')Gi(x, x'; z). (3.19) 

This integral kernel generates an bounded operator (use the estimate p.ip and employ the 
Schur- Holmgren criterion). We denote this operator by Sb{z). There are two important facts 
related to this kernel. First, Gi(x, x'; z) solves the distributional equation 

(P2 (0) + l^(x) - z)Gi(x, x'; z) = 5(x - x'). 

Second, one can prove that Sb{z) maps Schwartz functions into Schwartz functions. Moreover, 
employing (|3.18p and integrating by parts, we can establish an identity which holds at first in the 
weak sense on Schwartz functions: 

{{H^{B) - z)*, Sb{z)-E) = (1 + B Tb(z))S), (3.20) 

where Tb{z) is the operator generated by the following integral kernel: 

Tb(x,x';z) e'-^*^" (3.21) 

• I - 2iA(x - x') • VxGi(x, x'; z) + B| A(x - x')pGi(x, x'; z)} . 
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Using (|3.ip . and the fact that |A(x — x')| < |x — x'|, we obtain the foUowing point-wise estimate 
true for aU x 7^ x': 

max{ |5b (x, x' ; z) I , ITb (x, x' ; z) I } < const——— • (r) . (3.22) 

|x X I 

Clearly, Tb{z) can be extended to a bounded operator. 

Now let us prove that (j3.20p also holds true in the strong sense. Because Hoo{B) is essentially 
self- adjoint on the set of Schwartz functions, (|3.20p can be extended to any S S and any 
e 'Dom{Hoc{B)). It means that the range of Sb{z) belongs to the domain of Hoq{B) and 

{HUB)-z)Sb{z) = 1 + BTb{z). (3.23) 

At this point we can establish the following identity, valid for all z of interest: 

(iJoo(B) - z)-^ = Sb{z) - B{Hoo{B) - z)-^Tb{z). (3.24) 

Denote by Tb{z) := T^{'z) the bounded operator generated by the kernel 



Tb(x, x';z) := Tb(x',x;z). 
Now replace z with z in (I3.24|) and then take the adjoint. We obtain: 

(H^iB) - z)-i = Sb{z) - BfB{z){H^{B) - z)-\ (3.25) 

Denote by i^i,B(x, x'; z) the integral kernel of {Hoo{B) — z)^^ . With (|3.25p as starting point, 
together with (|5.5p . we can use the same argument as in the zero magnetic field case, in order to 
show that the resolvent (with exponential weights) maps into L°" . Thus its kernel obeys an 
estimate like in p.l5p . Moreover, (|3.25p implies via the Cauchy-Schwarz inequality that 

^|x— x'l 

|Xi,b(x,x';z)| < const ^--^^^ — • {r)"' (3.26) 



Moreover, we can repeat the arguments from the proof of Proposition l3.1l and get the boundedness 
and joint continuity for the kernel of Pj{B){Hoc{B) — z)^^ for large N, since all we have to do is 
to change Dj with Pj{B) and to notice that the formal commutator [Pj{B), HooiB)] is only linear 
in PkiBYs, therefore Pj{B) "commutes well" with {Hoo{B) ~ z)"^. 

Now let us show that {Hoq{B) — z)^^ maps into Holder continuous functions. In fact, one 
can prove the following estimate: 

Lemma 3.3. Fix a compact set U C M^, and fix j3 G (0,1/2). Take ijj with HV'Hl^ — 1. Then 
there exist two positive constants C and M such that 

sup {ry\{{HUB) - z)-V}(x) - {{HUB) - z)-'^}{y)\ 

<C.|x-y|^ (3.27) 

for any x, y G U . The same estimate holds true for Sb{z). 

Proof. The domain of Hoo{B) is locally , hence for some positive A, the function (iJoo(S) + 
A)~^'0 is locally . From the Sobolev embedding lemma, we obtain that {Hoo{B) + A)^^?/" is 
/3-II61der continuous for every /3 e [0, 1/2). The estimate (|3.27l) follows from the resolvent identity 

(iJoo(B) - z)-' = {Hoo{B) + A)-i + (z + X){HUB) + Xy^HooiB) - z)-\ 

The same result for Sb{z) follows from (|3.24p . □ 
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We now can start the actual proof of (ii) . We integrate by parts N times with respect to z in 
the expression of Fao, and N is supposed to be large. The terms we obtain in the integrand will 
look like this one: 

P,{B)iHoo{B) - ziy'Pk{B)iHoo{B) - Z2)-'^\ (3.28) 

where Ni + N2 = N + 2, and either A^i or N2 is large. Here zi and Z2 are complex numbers like 
z G or z ± Lo. 

By repeated commutations, we can always write this operator as {Hoo{B) — zi)~^W {Hryo{B) — 
22)"^, where is a sum of terms like this one: 

Wi := {H^{B) - zi)-'''v^,P^,{B){H^{B) - z^)-^-v^,P^,{B){Hoo(B) - z^^K 

Here ni,n2,n3 > 1, fa's are smooth and uniformly bounded functions. Note that such a term 
always starts and ends with a resolvent. Using the fact that Pai{B){Hoo{B) — zi)~^^Pa2{B) is 
a bounded operator, then we can always write W as the product of the form {Hao{B) — zi)^^W 
where is a bounded operator whose norm increases at most polynomially in r. Thus e^^''^^^'^W 
is a Hilbert-Schmidt operator, with a H-S norm which increases polynomially in r. Hence it is an 
integral operator which obeys the estimate: 

-^(y)^/^\W{y,y';z)\^dy dy' < const r*^ (3.29) 

where the constant and M are independent of r. Thus we have: 

Pj{B){H^{B) - z,r^'Pk{B){H^{B) - Z2)-^^ 

= {H^iB) - zi)-ie<->*>e-<->*>M^(i?oo(S) - ^2)"'. (3.30) 
Using an identity like the one in (|5.8[) . we can rewrite the above operator as: 

e<->*/'^(i7oo(S) - z^r'Te-^-^'^'-WiH^iB) ~ z^y^ 

where T is a bounded operator uniformly in r if (5 is small enough. Thus W := Te~^'^*/'"W^ is 
Hilbert-Schmidt, and has an integral kernel whose norm in L^(R^) is polynomially bounded in r. 

Therefore we are left with the investigation of the joint continuity in x and x' of the integral 
kernel defined by: 

/oo(x,x') := J y"ii'i,B(x,y;zi)W^'(y,y')ifi,s(y',x';22)dydy', 

where J J \ W {y,y')\'^dy dy' < const r*^. Using this, (|3.26p . and the Cauchy-Schwarz inequality, 
we obtain |/(x,x')| < const "r^^^ ||M^'||b2j uniformly in x,x' G U. Since W is Hilbert-Schmidt, we 
can approximate W^'(y,y') with a finite sum of the type ^gj{y)hj{y') where g's and /I's are L^. 
Thus we can (uniformly in x, x' G U) approximate the function /oo with functions of the type 

n 

Y,{{H^{B) - zi)-i.9,}(x){(7?oo(S)-Z2)-%}(x'), 



which from Lemma 13.31 we know they are continuous on compacts. Hence /oo is jointly continuous 
on its variables. As for the integral kernel of i^oo, we see that it can be written as the integral 
with respect to z of a finite number of kernels of the same type as /oo- Due to the exponential 
decay of fpD (see (jl.7l) ). wc see that we can approximate _F'oo(x, x') uniformly on compacts with 
continuous functions, and we are done. 

Therefore the function sb as defined in Theorem 1 1.11 is a continuous function. If V is periodic 
with respect to Z'^, then Hoo{B) commutes with the magnetic translations, defined for every 
7 e as (see also (IXTS)) ): 

[M^V](x) := e*-^'^o(^''^V(x - 7), Af-^M^ = 1. 
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Hence we have that as operators, M^^FooM^ = Fqo, which for kernels gives 

g-*B<^o(x,7)jr(x + 7, x' + 7)e'-^'^''('''+''''')j^(x, x'). (3.31) 

Now since (^0(7,7) = 0, when we put x = x' we get sb(x + 7) = sb(x) and we are done with 
(ii). □ 

4 Proof of iii. 

We start by proving the thermodynamic hmit for the conductivity, that is (|1.12|) . We need to 
introduce a special partition of unity in A^. 

4.1 Partitions and cut-offs 

Fix < a < 1 (small enough, to be chosen later) and define for t > 0: 

EL{t) := {x e Al : dist{x, OAl} ^tL"}. (4.1) 

This models a "thin" compact subset of A^, near the boundary, with a volume of order Be- 
cause we assumed that the boundary 9Ai was smooth, all points of Si(t) have unique projections 
on 9Al, if L is large enough. 

Then if ti < t2 we have Si(ti) C ^£(^2) and : 

dist{SL(il),AL\SL(i2)} > it2 - <i)i". (4.2) 

The subset Ai\Si(l) models the "bulk region" of A^, which is still "far-away" from the 
boundary. 

Now consider the inclusion of Sl(2) in the dilated lattice L"'!? . That is we cover Sl(2) with 
disjoint closed cubic boxes parallel to the coordinate axis, centered at points in L"Z^, of side 
length L". Denote by i? C L"^!? the set of centers of those cubes which have common points with 
Sl(2). Clearly, due to volume considerations, #i? ^ L^^^". 

In order to fix notation, let us denote by -^^(7, s) the cube centered at 7 e E', with side length 
equal to s > L". Moreover, denote hy E :— Eyj {(0, 0, 0)} (note that the origin cannot belong to 
E \i L \s large enough). 

Now choose a partition of unity {g-y\^^^§; of A^ which has the following properties: 



supp(.go) C Al\^l(1); (4.3) 

supp(g^)cX(7,2L"),7ei?; (4.4) 

0<5^<1, ^5^(x) = l, VxeAi; (4.5) 

||C''57lloo~i-"l^l, V/3eN3, y^eE. (4.6) 



This partition has the property that if we restrict ourselves to E^ then uniformly in L, the number 
of (7-y's which are not zero at the same time is bounded by a constant. Only 50 has ~ L^~^" 
neighbors whose supports have common points with supp(5o)- 

Now we choose another set of functions, {g-f^-^^E having the following properties: 



supp(go) C Ai\Si(l/4); go(x) = 1 if x e A^ \ Sl(1/2); (4.7) 

supp(g^) c X(7, lOL"), 7 e i?; g^(x) = 1 if X e K{^, 9i") (4.8) 

1 11?% I loo --i"""'', V/3eN^ y-fC^E. (4.9) 
These functions were chosen "wider" than the g-y's, and obey 

9i9-t=9i^ dist{supp(L'g^),supp(g^)} - L", 7 g i?. (4.10) 
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By Dijj we mean that we take at least one derivative of g^. 
Now let us define the third type of cut-offs, {g^yj^E- 

suM9j)cK{j,12L"),jeE; 5^(x) = 1 if x € i^(7, Hi") (4.11) 

Ili'^I^IU'-i"""'', V/5GN^ y^eE. (4.12) 

Note that we have g^g-y — g-y (the origin is not considered here). 

4.2 Proof of f fTTT^D 

Define for every 7 e -E: 

A^(x) :=|^(x) A(x-7). (4.13) 
Due to the support properties of our cut-off functions, we have the estimates 

57 =57 A(x-7), 

l|A^llci(R3) < const • L". (4-14) 

Define for every j G E (see also (|1.3p and (|1.4p ): 

P^{B) :=P(0) + BA^, i?L(S,7) ■.^P^{B)^ + V, (4.15) 

where the Hamiltonian is defined with Dirichlct boundary conditions. Note that Hl (i?, 7) — Hl (0) 
is a relatively bounded perturbation of Hl{0). 
Define the operator 

Ul{B,z) ■.^goiH^iB)~z)-'go 

+ ^ e^BM-n)g^(^HL{B,j) ~ z)-ie-^^*°(-'^)57- (4-16) 

One can prove that the range of Ul{B, z) is in the domain of Hl{B) and we have: 

{Hl{B)-z)Ul{B,z) = 1 + Vl{B,z), (4.17) 

Vl{B,z) {-2*(V5o) •P(B) - (Ago)}(ffcc (S) - ^)-\go 

+ ^ e^BM-n){_2^(y~g^) . V,{B) - {A~g,)}{HUB,j) - zr^e"^^'^''^-''^ g,. 

In order to obtain this equality we used the locality of our operators and various support properties 
of our cut-off functions, the identity (|3.18p . definition (|4.13l) . and (|4.5p . 
Then we can write 

(HLiB) - z)-' ^ Ul{B, z) + (HLiB) ~ zy^LiB, z). (4.18) 

The good thing about Vl{B,z) is that its operator norm is exponentially small in L". This is 
because we have the boundedness from (|5.5p and (|5.6p (valid also for Hl{B,'^), as can easily be 
seen from the proofs), and because of the estimate in (|4.10p . Indeed, for terms involving 7 7^ 0, 
put Xq = 7 in the two exponential estimates, and take +d on the left and —S on the right. Then 
we gain an overall decaying term from the left as (|4.10p implies: 

sup sup e-T^(<""^>"<"'"^>' < e-^^°, 7^0, (4.19) 

xGsupp(Dg^) x'Gsupp(gT) 

where Si > is small enough and L is larger than some Lq- 
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For 7 = the situation is slightly different, because we did not assume convexity for A^. But 
one of the terms whose norm we need to estimate is (see (|4.17p ) 



{W~go)-P{B){H^{B)-z)-^go. 
From p.ip foUows that the integral kernel of this operator is bounded by 

|(V5o) • P{B)iH^.{B) - z)-igo|(x,x') < const(r;)(r)^^e-T^I''-"'l. (4.20) 

Because x and x' are always separated by ~ (see the support properties for our cut-offs), we 
can write: 

l(Vgo) • P{B)iH^{B) - z)-i.9o|(x,x') < const(,7)(r)*^e-^l^-'''le-^^°, (4.21) 

where 5i and 62 are smaller than 6. 

Therefore we can write for all > 1: 

WVUB, z)\\ < const(77) • i^-Za^^^A/g-^L" ^ const(r?, a, N) ■ L-^{r)^^' , (4.22) 

where we have to remember that we have ~ L^~^" of g-y's. The second estimate says that the 
norm decays faster than any power of L, with the price of a higher power in (r). 

We now want to show that Vl{B, z) does not contribute to the thermodynamic limit of aL{B). 
We have the following result: 

Proposition 4.1. 

aL{B) = - ] -Tr / fFD{z){Pi{B)UL{B,z)P2{B)UL{B,z + Lo) 
VoI(Al) Jt^ 

+Pi{B)Ul{B, z - uj)P2{B)Ul{B, z)} dz + 0{L-°°). (4.23) 

Proof. The main idea is to show that when we replace {Hl{B) — z)^^ in Fl with the right hand 
side of (|4.18p . all terms containing Vl{B) will generate (after integrating by parts with respect to 
z) some operators which in the trace norm will decay faster than any power of L. 

If the differentiation does not act on Vl{B) but on the other resolvents, then it is enough to 
know that in the norm of B{L^) it goes to zero faster than any power in L, as we saw in (|4.22p . 

If the differentiation with respect to z acts on Vl{B), then there will be a few terms which 
must be separately considered, and prove their smallness in the trace norm. 

To give an example, after differentiating A^ — 1 times with respect to z (N large) , we obtain a 
term containing a factor like: 

3 

Y,idf9o)Pj{B)iH^iB) - z)-%. (4.24) 
i=i 

We will prove (in the trace norm) that it decays faster than any power of L, times some polyno- 
mially bounded factor in (r). Note that all the other factors multiplying the above operator are 
bounded operators, with a norm which is polynomially bounded in (r). 
Let us start with a technical result: 

Lemma 4.2. Let Qi and Q2 be two compact unit cubes such that dist((5i,Q2) — d > 1, and let 
Xi, Xi denote their characteristic functions. Let a S {0,1} and j € {1,2,3}. Then if N is large 
enough, there exist three constants S2 > 0, Ni > 1 and C > 0, all three independent of z G T^^, d, 
a, j and Q 's such that 

\\XiP^{B){H^{B) - z)-^X2||bi < Cr^^ exp{-d<52/r}. (4.25) 
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Proof. We assume that a — I, the other case being similar. The strategy is to write our operator 
as a product of two Hilbert-Schmidt operators. By commuting Pj{B) with one resolvent, we can 
rewrite our operator as: 

XiiHooiB) - z)-'T,iHoo{B) - z)-^+\2, 

where Tj is a bounded operator which contains factors like Pk{B){HoQ{B) — z)^^. 

Denote by X2 an arbitrary point in the support of X2- We insert some exponentials in the 
following way: 

li $2 < 6 we can write: 

g-^(x-X2)g^(x-x')g^(x'-X2) < ^^^^^g-l^|x-x'| 

Now the two factors containing resolvents in (|4.26p are Hilbert-Schmidt due to the presence of 
the cut-offs X and the exponential decay of our kernels (Proposition 15 . 21 and (|3.ip ). Their Hilbert- 
Schmidt norm will grow polynomially with r, but be independent of d. The choice of X2 provides 
the decaying exponential factor on the right hand side of (|4.25l) . □ 

Let us go back to (|4.24p and try to use the previous lemma. We will show that in the trace 
norm, this operator decays exponentially in L". Consider only one j. Cover both supp{i9jgo} 
and suppjgo} with disjoint cubes centered at points in Z'^ and side length equal to 1; we only use 
cubes which have common points with the respective supports. We then have 

{d,~go)Pj{B){H^{B)-z)-''g„ 

= Y,Xs{dj~9o)Pj{B){H^{B) - z)-^goXs', (4.27) 

where Xs and respectively Xs denote the characteristic function of such unit cubes which cover 
supp{djgo} and respectively supp{.go}. The number of cubes needed to cover the support of go 
is of order L^, while for the other one is of order L'^"; hence we have about L^+Sa ^gj-j^g the 
above double sum. 

But each operator of the form 

Udj~9o)PAB)iHociB) - z)-%Xs' 

is exponentially small in the trace norm due to Lemma 14.21 since the distance between any two 
supports of Xs and Xs' is of order L". Hence the entire sum in (|4.27p will be (r dependent) expo- 
nentially small in L". But then we can trade off the fading exponential decay with a polynomial 
decay in L and a polynomial growth in r as we did in (j4.22p . So this term is under control. 

Now let us go back to the beginning of the proof of Proposition l4.1l Other "bad" terms from the 
remainder in (|4.23|) after differentiation with respect to z will contain powers of {Hl{B, 7) — 2)^^, 
like for example 

{dfg^)Pj^^{B){HL{B,j) - z)-^g,. (4.28) 

Here we cannot easily commute with P's due to various boundary terms. But we do not need 
to do that. Look at (|5.8|) . where we put s — and Xq = 7. Because is bounded from above 
by (see (jiTi]) ). it follows: 

IK-Az, + l)(i/L(S,7) - < {r)''L^- (4.29) 

It means that for small S, the resolvent e''^^'~''\HL{B,j) — z)~'^e~''^^"^'''^ sandwiched with 
exponentials remains Hilbert-Schmidt (with a norm which does not grow faster than the B2 norm 
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of {-Ad + 1)"^ times and some polynomial in (r)). Since we have N resolvents, the product 
of two of them will give a trace class operator. Now we can repeat the insertion of exponentials 
as we did in (|4.26p . and use (|4.19p for getting the exponential decay in L". 

□ 

Now let us show that all terms involving the sum over 7 G -E in (|4.16p will not contribute at 
the end. The explanation is that these terms are "localized near boundary". We can formulate 
the result as follows: 

Proposition 4.3. If a > is small enough, then: 
lim \aL{B) + 

— i— • Tr / fpDiz) \Pi{B)UHoo{B) - z)-^gMB)UHoo{B) - z - uj)-^g^ 
VoI(Al) L 

+PiiB)go{HooiB) - z + Lor^goP2{B)go{H^XB) - z)-'gQ ] dz} = 0. (4.30) 

Proof. If we compare this with (|4.16p . we see that we need to show that all terms containing 
factors localized near boundary will converge to zero. Let us look at one such term, and prove the 
next lemma: 

Lemma 4.4. Assume that < a < 1/3. Then we have that 

voik) ■ ^'L ^--(^)^^(^)^°(^-(^) - ^)"^° 

• P2iB) J2 e'^'*"'^-^^^~g^{HL(,B,-i) -z- uiY^e-'^'*"'^-^^^ g^dz = 0. (4.31) 

Proof. Using (|3.18p . and the fact that 5^57 — g-y, we can rewrite the above term as 
TT^, • E / fFD{z)Pi{B)UHoo{B) - zY^g^l^ 
■ e'^'^°^-^^^P2^^{Byg^{HL{B,-i) -z- uY^e-"^'^''^-^^^ g^dz. (4.32) 



After integrating by parts iV — 1 times with respect to z, we have to deal with several situations. 
Let us take one resulting term (just the operator in the integrand): 

Pi{B)UHoo{B) - z)-i5ol^e*^*«(-'^)P2,^(S)ff^(i?L(S,7) - z)-^e-^^^°(-''^)<?^, (4.33) 

and let us estimate its trace norm. The factor containing Hoo{B) is just bounded, and we cannot 
use it as a Hilbert-Schmidt factor. What we do is to commute g^ over one resolvent to the right 
and get the identity: 

P2,^{B)g^(HL{Bn) - z - u^g^ (4.34) 
= P2^^{B){Hl{B, 7) - 2 - Lo)-'g^{HL{B, 7) - z)-^+'g^ 
+ P2,^{B){Hl{B, i)-z- u;)-\Hl{B, 7), g^]{HL{B, - z - co)-%. 

Now the second term will again contain at least one derivative of g-y and reasoning as we did for 
(I4.28P we can show it will be exponentially small. Let us look at the first term. The operator 
P2,-i{B){H l{B ,"f) — is only bounded. But g^{HL{B,^) — ()~^^^g'y is trace class and 

\\~g,{HUB,j)-Cr''+'9,\\B, (4.35) 
< -l-^\\~g^{HL{B,^) - C)-'\\b. ■ \\{Hl{B,^) - Cr'a-rllB.- 
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Hence using (|4.29p we have 

\\g-,{HL{B,^) - Cr^Ws, < const(r)*^||g^(-Az5 + 1)-'\\b, ■ i". 

But the Hilbert-Schmidt norm of g^{—A]j + 1)^^ is of order of the square root of the support of 
g^, that is L-^"/^ (use here (|5.13p with B = 0). The other factor comes with a similar contribution, 
hence we can write 

Wg^iHUB,^) - z - c.)-^+\g^||B, < const (r)'^ ■ L^^ . 

Now using this in the integral with respect to z, this particular term will give a contribution of 
j^5a Jqj. each 7. Since we have ~ different 7's, the total contribution will be bounded by 

jji+^a ^ But if a < 1/3, after we divide with the volume of it will converge to zero. 

Now let us go back to (I4.32[) , and see that after integration by parts we can get a term like 

Pi{B)UHoo{B) - z)-"" g^l^ 

. e^BM;^)p^^^^Brg,{HL{B, j) - z ^ u;)-'g,. (4.36) 

Here the operator P2^^{B)g^(HL{B,'j) — z)^^ is not Hilbert-Schmidt, so we have to look at the 
first factor. 

We can write 

Pi{B)go{H^{B)- z)-^go~g-, (4.37) 
= Pi{B)go{Hoo{B) - zy+%{HUB) - z^'go 
+ Pi(B).go(i?oo(S) - z)-''[H^iB),~g^]{H^.{B) - z)-^go. 

In the first term, the function g^ makes the two resolvents next to it become Hilbert-Schmidt, 
each having a norm proportional with L^"/^ and some power of (r) (use the exponential decay of 
the kernels). So this term is "good", considering that we have to divide with in the end. The 
second term contains at least one derivative of g^ (here [Hao{B), g^] is linear in Pj{B)'s), together 
with factors like 

{P^{B)UHoo{B) - z^ P,(B)}{dfg^)[H^{B) - z)-'go. 



Now here we can use the estimate from Proposition 13.21 and see that we again have a product 
of two Hilbert-Schmidt operators: the first Hilbert-Schmidt norm will be proportional with L'^/^, 
while the other one will behave like L'^"/^. When we divide by L^, this contribution will go to 
zero. 

All the other terms resulting from integrating by parts with respect to z can be treated in a 
similar way. □ 

Now let us go back to (|4.23p and analyze another boundary term: 
Lemma 4.5. For every < a < 1/3 we have: 

L^oo VoI(Al) Jr^ 

■P,{B) J2 e''''^°^-''''^~gy{HL{B,j') - z)-'e-'''*°^-''''^gy 

■ P2{B) e'-^'^°(-'^)5^(i?L(B, 7) - z - uj)-^e-'^''">^-^^'>g^dz = 0. (4.38) 

Proof. Let us note that when keeping 7 fixed, only a finite number (L- independent) of 7"s will 
have an overlapping support. This means that the above double sum will only contain around 
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lJ2-2a i^on-zero terms. Now use again (13.17^ and integration by parts with respect to z. Each 
non-zero term in the double sum wiU be a product of two Hilbert-Schmidt operators, each with a 
Hilbert-Schmidt norm of the order of L^"/^. The total trace norm will grow at most like L^+'^", 
hence if a < 1/3 this term will not contribute. We do not give more details. □ 

The last ingredient in proving (|1.12[) is contained in the following result (see also (14.30^ ): 
Proposition 4.6. 

1™ -^^TTTTT • / fFDiz) {Pi{B)go{H^{B) - zr^goP2{B)go{HUB) - z - Lo)-^go 
L^oo VoI(Al) Jy^ ^ 

Pi{B)g~oiH^{B) - z + Lo)-'goP2{B)g~o{H^iB) ~ z)-^go} ^z ^ - J^^ SB(x)dx. (4.39) 

Proof. First, due to support properties, we have 

goP2{B)g~o = gMB) = P2{B) - (1 - go)P2{B). (4.40) 

Because 1 — go is supported outside of a thin region around the boundary of A^, , all terms generated 
by (1 — gQ)P2{B) will converge to zero; let us prove this. After integrating — 1 by parts with 
respect to we obtain several terms from the integrand which look like this one: (A^i + N2 > 
N > 10): 

Pi{B)g~o{Hoo{B) - z)-^Hl - go)P2{B)iHoo{B) - z - coy- go. (4.41) 

Due to the symmetry of this term, assume without loss of generality that A^i > 5. Then by writing 
1 — go = (1 ^ .9o)XAi + (1 — .9o)(l ~ XAl)i get two types of contributions. The one coming from 
(1 — go)(l ~ XAl) ~ (1 ^ XAl) is localized outside A/,, and its trace norm will be exponentially 
small in L"; we can apply Lemma 14.21 since the distance between A'^ and the supports of go and 
go is of order L". The number of needed covering unit cubes is polynomially bounded in L. 
Thus the only contribution from (|4.4ip can come from: 

Pi{B)g~oiH^{B) - z)-^i(l - go)xAMB){HUB) - z - uj^-go- (4.42) 

Here we can write: 

Pi{B)go{Hoo{B) - zr'^'il - .go)XA. (4.43) 
{P^{B)go{Hoo{B) - z)-^} ■ {{HUB) ~ z)-^^+2(l - go)XAj, 

where both factors are Hilbert-Schmidt, with kernels exponentially localized near diagonal as in 
Propositions 13.11 and 13.21 The Hilbert-Schmidt norm of the first factor is bounded by i^/^, while 
for the second one we have a bound of L^+^Z^ (square roots of certain volumes). The trace norm 
of the product is thus bounded by L^/^+^Z^ and some polynomial in (r). After integrating with 
respect to z, and dividing with (the volume of A^), this term will converge to zero provided 
a < 1. 

Now we can go back to (|4.40l) and analyze the term generated by P2{B). Because there are no 
other cut-offs in the middle, and because the commutator [Pi{B),go\ will generate another fast 
decaying term, we see that we have just proved the following identity (see Theorem 1 1.1 1 ii for the 
definition of Fao)'- 

ii-^ (^^(^) + Volk)^^^^°^-^°0 ■ (^-^^^ 

But the operator goFoogo is trace class, with a jointly continuous kernel, hence (see (|4.ip . (|4.3p 
and (03)) 

Trjgo^'oogo} = / SB(x,x)go(x)dx. 

Jsupp(go) 
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Using the periodicity of sb with respect to 1/^ and the support properties of go, we finally get: 

lim <j]^{B)^~ / SB(x,x)dx (4.45) 

and the proof of (|1.12p is over. 

□ 

4.3 Proof of dlllSl). 

We start by investigating dso'LiO) and try to put it in a form which is better suited for the 
thermodynamic limit. 

Note that we have already argued that cr^ was smooth in B near zero (see the remark around 
(|2.5p ). The hard part is to show that the polynomial growth in L of the trace norm does not 
appear in the actual trace. Similar difficulties involving magnetic semi-groups were encountered 

in miiniiiiT]. 

4.3.1 Only Ul{B) counts. 

Notation. In order to shorten our formulas, we will sometimes write {z ^ z — lo), which means 
that we have to repeat the previously appearing product of operators in which we have to change 
z with z — u. 

First, let us prove that even if we differentiate with respect to B, we still have a result similar 
to Proposition 14. II 



Proposition 4.7. At B = Q we have: 

Ob |fTL(B) + yj^^^^ ■ Tr^ fFoiz) {Pi{B)Ul{B, z)P2(B)Ul{B, z + oj) 

+ z^ z-uj}dz^ =0{L-°°). (4.46) 

Proof. As usual, before differentiating with respect to B one has to perform a certain number of 
integrations by parts with respect to z. There will appear many terms in the remainder containing 
Vl{B, z) (see (|4.17p V but all of them will have the distinctive feature of containing pairs of cut-off 
functions whose supports are at a distance ~ one from another. 

We do not want to treat all the situations which can appear here, and instead we will only 
prove a typical estimate needed for the result. 

Lemma 4.8. Assume that N is large enough, and choose z such that dist{2:, [0, oo)} — r/ > 0, 
r = 'Si{z). Then the map 

(-1, 1)3B~. (A5o)(i?oo(B) - z^go G Bi{L^) 

is differentiable in the trace norm, and there exists di small enough and M large enough such that 

\\{dB{Ago){H^{B) - z)-%}b=o||i ^ const(77, 7V)e-^^° (r)*^ (4.47) 

Moreover, 

\\{Ago){H^{B) - z)-^go - {A~go){H^{0) - z^ go 

-B{9B(Ago)(i?oo(S)-z)-%}B=o|L, < const(r;, 7V)e-^^° (r)*^ (4.48) 
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Proof of the lemma. We only prove differentiability at i? = 0. The key ingredient is to give a 
proper sense to the formal identity: 



= -{H^{B) - z)-^{Hoo{B) - i/oo(0)}(-ffoc(0) - z)-\ 



(4.49) 



As it stands, the right hand side makes no sense because Hoo{B) — H^{0) contains terms like 
— 2iBA(x) • Vx which are not relatively bounded to the free Laplacian due to the linear growth 
of our magnetic potential. For < 82 < 5^ and Xq G M^, we can write (still formally): 



(4.50) 



-eT^(-^°\H^{B) ~ zy^{H^{B) - H^{Q)}{H^{Q) - z^^e 



Note the fact that the growing exponential is weaker that the decaying one. This identity now 
holds in the Hilbert-Schmidt class; in order to see this, introduce the notation 



82' 



(4.51) 



Then we can write the above identity as 

[e^^-'^^'^H^iB] - i/oo(0)}e-^<-''">} 
>(iJ^(0)-z)-ie"^<---)} 



(°l-l)^2 , 



(7) (--xo; 



(4.52) 



Now Hoo[B) — i?oo(0) = 2i3A • P(0) + B^A^ is proportional with B, and the linear growth of A 
is compensated by the higher exponential decay on the right hand side. Hence (15. 5|) and ()5.9p (at 
B = Q) imply that if 82 is small enough, we have: 



^^■"""^{i/oo(i?) - i/oo(0)}(ffoo(0) - z)-ie-^<-''°> 



< const (xq) 



\B\ ■ (r 



M 



(4.53) 



The estimate p.26p (valid uniformly in _B e [—1, 1] and the S there should be chosen slightly larger 
than the S3 ) tells us that: 



Hence we have proved the estimate 



< const (r) 



M 



;^<-"°>{(i?oo(B) - z)"i - (ffoo(O) - z)-i}e- <*<-^''°> 

Now go back to (|4.50p and isolate the linear term in B: 

e^(-^"\{H^{B) - z)-i - (i/oo(0) - zr'je-^'^-^"^ 
= -i3e^<— o)(i/^,(0) - z)-'{2A ■ P(0)}(ffeo(0) - z)-' 
- B^e^''-^"\H^iO) - z)-i{A2}(i/^(0) - z)-'e-^ 



< const|_B|(xo 



|2 (r)^^. 



(4.54) 



(4.55) 



(4.56) 



(-xo> 



eTA<— o){(if^(i3) - z)-i - {H^{0) - 2)-i}e-T^<-^"> 



■ e <-■> 



'^-'"'^ {H^{B) - ffoo(0)}(ffoo(0) - z)-ie-T*<-''«\ 
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where in the last term we inserted an exponential with 62 < S < 63. Now it is easy to get the 
estimate: 

+Be^<-"°>(i/oo(0) - z)-'{2A ■ P(0)}(i?oo(0) - z)-ie-^<-'"''> 

B2 

= const|B|2(xo)MO*'^- (4.57) 



This is enough to prove that if iV > 2, the mapping in Lemma [4~8l is differentiable in the trace 
norm sense at B = 0. Indeed, proceeding as we did for (|4.27p . we can insert many cut-off functions 
and cover the supports of go and go. Using the same notations, we have that for some S2 < 63 < Si 
we have 



Xs{^~go){H^{B) - z)-'' goxs' (4.58) 



Now we can differentiate with respect to B in the trace norm-sense, using the result for the 
Hilbert-Schmidt norm and the fact that we have at least two factors in B2 (adapt (|4.54p '). At last 
we again use the fact that the distance between the supports of x's is ~ L", and that all growing 
factors are just polynomials in L. We consider Lemma 14.81 proved. □ 

Now we can use (|4.48p in all the terms on the right hand side of (|4.46p which contain operators 
of the type treated in Lemma 14.81 The exponential decay of fp d can be used to obtain a decay 
faster than any power of L. 

All other terms from the remainder can be treated in a similar manner, and we consider 
Proposition [4?7] as proved. □ 



In the remaining part of our paper we will investigate the thermodynamic limit of the main 
contribution to dsO'LiB), given by 

9b Tr^ fFDiz){Pi{B)ULiB.,z)P2iB)ULiB,z + Lo) + iz^ z-Lu)}dz. (4.59) 
4.3.2 The boundary terms 

Here the magnetic perturbation theory will play a crucial role. There are several terms in the 
definition of Ul (see (|4.16p ). and when we insert them into (|4.59p they will generate even more 
terms. 

We will now prove that only the term which contains two resolvents with Hoc will contribute 
at the thermodynamic limit. The main result can be stated in the following way: 

Proposition 4.9. At B = and a sufficiently small we have: 

■ {Pi{B)UHoo{B) - zr'goP2iB)go{H^iB) - z - uj)-'go 

+ z ^ z - I = 0. (4.60) 

Proof. Let us start with the following boundary term: 
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Xi{B,L) :-Tr^ fFD{z) { PiiBygoiH^iB) - z)-' go 
+ (z ^ z — tj) I dz. 



9i 



(4.61) 



We have already seen in Lemma [4.41 that Xi{B,L) divided by the volume of converges to 
zero when L converges to infinity. Now we would like to show that {9bXi}(0, L) has the same 
property. 

Proceeding as in (|4.32p . we can rewrite Xi{B, L) as 



Xi{B,L,j) := Tr^ fpoiz) { Pi{B)go{HUB) - z)-^go 

+ {z — > z — cj) I dz. 

We now prove the following estimate: 
Lemma 4.10. For every L > 1 we have 

1 



(4.62) 



(4.63) 



sup sup 

7e-E0<|B|<l 



B 



{Xi(B,L,7)-Xi(0,L,7)} 



< const(a) • L 



5q 



(4.64) 



Remark. Before starting the proof of this lemma, let us note that it immediately implies that 

1 



VoI(Al) 



{dBXi}{0, L) ^ 0. 



(4.65) 



This is so because we know that {(?bXi}(0, L) exists, and moreover, it must be bounded from 
above by the right hand side of ()4.64p times the number of 7's in E, i.e. ~ Then if a is 

chosen small enough, after dividing by ~ we get something converging to zero. 
Proof of Lemma I4.10i Define the function 

1i(B,L,7) := Tr^ fpoiz) { Pi{B)goSB{z)go 
■ e^-^'^°(-''')F2(0).g-,(i?L(0) -z- u;)-^e-'^'^"'^--^'^g^ 

+ {z^ z-Lu)^dz. (4.66) 

The proof of this lemma has two parts. The first one will state that 

1 



sup sup 

7GE 0<|B|<1 



B 



{Xi(i?,L,7)-Xi(B,L,7)} 



< const(a) • L 



5a 



while the second one is 



sup sup 

7e-E0<|B|<l 



-{Xi(i3,L,7)-Xi(0,L,7)} 



< const(a) • L 



5a 



(4.67) 



(4.68) 



Part one. The first estimate is not very much different from what we have already done until now. 
One uses repeated integration by parts with respect to z in Xi{B, L, 7) and then we expand each 
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resolvent {Hoo{B) — ()~^ using the first equality in (|3.24p . and each resolvent {Hl{B, 7) — C)"^ as 
perturbation of {Hl{0) — ()~^. Note that due to (I4.14p . the growth induced by A^(x) does not 
exceed L". As for the Hao{B), commuting Pj{B) with the magnetic phases will always transform 
A(x) into A(x — x'), whose growth will now be tempered by the exponential decay from p.22p . 
The general strategy for all terms arising from integration by parts with respect to z is to estimate 
the trace by the trace norm, using the fact that the trace norm of a product is bounded by the 
Hilbert-Schmidt norm of two factors. Several other estimates are needed in order to prove (|4.67|) . 
Remember that GAr(x, x'; z) is the integral kernel of (-ffoo(O) — z)~'^ , N > 1. Denote by s'^j^^z) 
the operator corresponding to the integral kernel e»-Bvo(x,x )(7^(x, x'; z). If = 1 they coincide 
with the operators Sb{z) defined in (l3.19|) .Then we have: 

i. For every Be [-1, 1], and > 1 

\\s\^\z)\\ < const (Ar)(?y^ (4.69) 

ii. For every B e [—1,1], N > 1 and Q C M'^ a compact set: 

IIXQ^lf ^(z)IIb. < const(iV)VVoI(Qy {r)'' . (4.70) 

iii. For every Be [—1, 1], 

\\{H^{B) - z)-i - Sb{z) + BSb{z)Tb{z)\\ < const|B|2 (r)^. (4.71) 

iv. For every B G [—1, 1], > 2 and Qi,2 C M.^ two compact sets: 

\\XQrS^B^\z)xQA\B, < const(7V)VVol(Qi) /MQ^ {r)^' . (4.72) 

The first and the third ones are easy consequences of (|3.24p . For the second and fourth ones we 
have to differentiate — 1 times in (|3.24p and write: 

4^)(z) = i-lf-\H^.iB) - z)-^ (4.73) 
-BJ2 hlf-'-'{Hoo{B) - z)-^+'=Ti'=)(z). 

Part two. We will now concentrate on (|4.68p . which needs a new idea. 

An heuristic argument. First we perform some formal computations, in order to illustrate how 
magnetic phases will transform the trace into a more regular object. Assume that the operator 
under the trace in (|4.66p has a jointly continuous integral kernel; remember that the operator 
Sb{z) defined in p.l9p had a magnetic phase. Commute this phase with Pi{B) as in p.lSp . and 
write the following formal expression for the integral kernel of the operator whose trace we want 
to estimate: 



dzfpDiz) \ / dx' 



g.i3^o(x,x')[p^^^(0) + i?Ai(x - x')]ffo(x)Gi(x,x', z).go(x') 
• e'^^°(^''^)F2,x,(0)ff7(x')(^L(0) - z - a;)-i(x',x")e-'^'^°(^"^^)57(x") 
+ (z->z-tj)}. (4.74) 

The above expression gives an integral kernel /(x, x"). If we could prove joint continuity, then we 
could write 

Ai(B,L,7)=/ /(x,x)dx. (4.75) 
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Now let us see what happens with the phases in (|4.74p when we put x = x". We have the identity: 

//(x,x',7) :-le3-[(x'-x)A(7-x')], 

</7o(x, x') + ^o(x', 7) - ^o(x, 7) + /^(x, x', 7). (4.76) 

The crucial thing is that when x = x" in (j4.74p . the magnetic phases cancel each other and only 
the flux fl remains. 

Thus the operator given by the following integral kernel 

^ dzfFD{z)S^ dx'e'^-^'('''''''^'[Fi,x(0) + SAi(x-x')]go(x)Gi(x,x',2)5o(x') 

• F2,x'(0)g7(x')(-ffL(0) - z - co)-i(x',x").9^(x") +{z^z-Lu)}, (4.77) 

must have the same trace since its kernel has the same diagonal value. Remember that this is just 
an heuristic argument, precise details are given in the next paragraph. 

The rigorous argument. We now start the rigorous proof of (|4.68p . We integrate by parts 
with respect to z in (|4.66p . and let us first focus on one term, namely the one obtained when all 
derivatives act on the resolvent in the middle: 



R{B,L,j) :=Tr /^^c(2) { Pi(i?)5o5B(z)go e*''^°(-^)P2(0)g-,(i?L(0) - z - c^)-^e-*^^«( '^)57 
+ {z^ z-u;)^dz, N>3. (4.78) 

Denote by {V-'fc}fc>i and {Xk}k>i the eigenfunctions and eigenvalues of i?i(0) respectively. 
Then the operator: 



•= E / fPDiz) I Pi{BrgoSB{z)go e*''^''(-'^)P2(0)57lV'fc)(^fcl7T — 

+ {z^ z-Lu)^dz (4.79) 
is trace class. Using 

Tr(i?K) - R{B,L,j)TrJ^ fMz) { Pi (5)30^5(^)50 e^^'^«(-^)P2(0).9^(iJi(0) - z - uj)-' 

■ E IV'fc)(^fcl71^ ^-^w^e-^^'^oi.Dg^ +(^ ^ ^ _ ^) I dz (4.80) 

we obtain 

lim Tr(PK) = P(P,i,7), (4.81) 

where we use the fact that the square of the resolvent is trace class (here — 1 > 2), together 
with the boundedness of the rest of the factors, with norms polynomially bounded in (r). 

Now let us show that Rk has an integral kernel i?i<-(x, x') which is jointly continuous on 
Al X A^. For, choose a point (xo,Xq) e x and let us prove continuity there. We know 
by elliptic regularity that tjjkS are smooth in A^. Then the function P2(0)g-yVfc is smooth in Al- 
Denote by xo a smoothed-out characteristic function of a small ball around Xg, whose support 
is included in A^. The operator goSB{z)go sends smooth functions into smooth functions, hence 
505*3(2)50X0^2 (0)57 "0*: is smooth in A^. Then since the integral kernel of Sb is smooth outside 
its diagonal, it means that gQSBiz)go{l — xo)P2iO)g'y'4'k is smooth at Xq, and remains so even after 
applying Pi(P). All bounds are growing at most like a polynomial in (r), hence the integral in z 
preserves the continuity. The variable Xq only meets smooth functions in A/,, and we are done. 
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We thus conclude that Rk is trace class and has a continuous integral kernel. For ev- 
ery increasing sequence of compacts fig such that /" (in the sense that il^ C and 
lim5_,oo Vol(Ai \ ils) = 0), we can write: 

Tr(i?K) = hm / i?K(x,x)dx. (4.82) 

Jo, 

Let us denote by Q-y{B, z) the operator given by the integral kernel 

Q^(x,x';S,z) (4.83) 

:= e*^/'(^'^'^^)[Fi,,(0) + i3Ai(x- x')]ffo(x)Gi(x,x', z)5o(x')^7(x'), 

At this point we can get rid of the magnetic phases using (|4.76p . and using the notation from 
(1^:551) we have 

Tt{Rk) - / fPDiz) I QyiB,z) P2{0)~g^\^p^){4,k\- 

~l JT^ V'^k — Z — Ul) 

+ {z^z-Lu)^dz. (4.84) 

Now Q-y is a nice bounded operator, and we can let ii" ^ 00, because the integral will converge in 
the trace class in the same way as in (|4.80p . Then (|4.81l) implies that: 

R{B, L,j)TtJ^ fFD{z) { Q^{B, z) P2{0)g^iHL{0) - z - c^)-% 

+{z^ z-ij)^dz. (4.85) 

Denote by Qy(B, z) the operator given by the integral kernel 

Q^(x,x';i?,z) := l(e^s/Hx,x',7) _ i)P, ,^(o)go(x)Gi (x, x', z)go(x' )^7(x') 

+ e'^^'(^'^''^Ui(x - x')go(x)Gi(x, x', z)go(x')^-r(x')- (4-86) 
We can write: 

R{B,L,^) - i?(0,L,7) = B Tr^ fMz) { QyP2{0)g^{HL{0) -z-uj^g^ 

+ {z^z-uj)y (4.87) 
We also note the estimates: 

|/;(x, x', 7)1 < i|x - x'l |x' - 7|, |e^s/'(-'-'.^) - 1| < B\fl{^, x', 7)|. (4.88) 

Now it is easy to see from ()4.86p . (|4.88p and (13. ip that belongs to B{L?'), with a norm bounded 
by (r)*^L". By writing 

{Hl{0) - z - Lu)-^g^ = (i7i(0) - z - uj)-^ ~g^{H L{i)) - z - uj)-^g^ (4.89) 
+ {Hl{Q) - z - u:)-\l - ~g-,){HL{0) - z - Lo)-'g^, 

we can see that the operator {Hl{Q) — z — uj)~^gy is trace class and 

\\{HLiO) -z- Lj)-^g^\\B, < const • {r)^'L^'', (4.90) 

after estimating the Hilbert-Schmidt norm of each factor in the two terms. The second one will 
be exponentially small due to the support properties of 5^'s. 
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We have thus proved 



sup sup 

7G-E 0<|S|<1 



i{i?(B,L,7)-i?(0,L,7)} 



< const • L"*". (4.91) 



After summation over 7, the bound is hke L^+^", and if a < 1/2 it wiU not contribute to the 
thermodynamic limit. 

Remember that this was just one possible term arising after integrating by parts — 1 times in 
(|4.66p . All other terms having sufficiently many derivatives acting on the resolvent, can be treated 
in a similar way. A different class of terms is represented by the one in which all derivatives act 
on iS's(z). Let us define: 

RiiB,L,j) 

:= Tr^ fFDiz) { Pi{B)~g^sf\z)go e'^^''(-^)P2(0)ff^(i?L(0) - z - L.)-ie-'^^°(-^)g^ 
+ {z-f z-uj)^dz. (4.92) 
We commute back ^2(0) over the phase at its left and write: 

Q2(S,^,7) := Pi{B)g^S^^\z)goP2{B)~g^, (4.93) 



Ri{B,L,^) = Tr fFDiz) { Q2(S,^,7) e*^^«(-'^)57(^i(0) " z - c.)-ie-^^'^«(-'^)g^ 
+ (z ^ z — oj) I dz. 

Now if A^ is large enough, by using (|3.18p . the regularity of Gn{x, x'; z), and the exponential decay 
of the kernels, one can prove an estimate (A^ large enough): 

|Q2(i?,2,7)(x,x')| < const(A^) • (r)*'^e-F) l''-'''l|^(x'). (4.94) 

It means that the integrand in Ri is a product of two Hilbert-Schmidt operators. We can again 
introduce the cut-off with the spectral projection of Hl{0), get rid of the magnetic phases and 
introduce the more regular phases, and so on. We consider that Lemma 14.101 is proved. □ 

Besides Xi treated in the previous lemma, there is only one other boundary term which needs 
special attention. This term is the one containing a double boundary sum: 

A2(B,L)= ^2(5,^,7,7'), (4.95) 

7,7'G-E 

A2(B,L,7,7') :=Tr^ fFuiz) 

. I e'BM:f')p^^^,^Brgy{HLiB,Y) - z)-'e-''"^o^-'''^ gy 
■ e^^'^«(-^)P2,7(i3)57(^fi(^,7) - z - c^)-ie-'^^«(-'^)57 

+ {z ^ z - uj)^ dz. (4.96) 

and we want to prove that for every i > 1 we have 

1 



sup sup 

7,7'e-B0<|B|<l 



^{A2(i?,L,7,7')-^2(0,L,7,7')} 



< const(a) • L"". (4.97) 



Note that this would again imply something like (|4.65p but for A2, because there is a finite, 
i-independent number of 7's and 7"s with joint support. 
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The strategy is the same. We define 

■ e'-«'^»(-''')F2(0)g^(i?L(0) - 2 - uj)-^ e-'^'^''^-^^^ g-y 

+ (z^z-cj) jdz, (4.98) 

and we want to prove two analogues of (|4.67p and (|4.68p . The analogue of (|4.67p is "easy" , 
but the analogue of (|4.68p again requires a limiting procedure which would allow us to write the 
X2{B, L, 7, 7') as the trace of a more regular object in B. The main point is that this new object 
will be given by the composition of those four magnetic phases present in X2. Namely, let us 
notice the following identity: 

(^o(x, 7') + ¥^0(7', x') + ¥'o(x', 7) + Ml, x) (4.99) 
= //(x,7',x') + /;(x',7,x) 

= //(x,7',7) + /;(x',7,7'). (4.100) 
Now (|4.100p allows us to write: 

X2{B,L,^,j') = Tr^ fMz) { e*^/'( '^''^)Pi(0)ffy(iJi(0) - z^gy 
■ e'^^'(-''^-^')P2(0)57(i?L(0) - z - LuY'g^ 

+ {z z -uj)^dz. (4.101) 

The good thing about this formula is that on the supports of 57's, these fluxes are at most of 
order L^", being bounded from above by |x — 7I • I7 — 7']. Remember that the non-zero terms 
must have I7 — 7'] < const • L". Now we can expand the exponentials and prove the analogue of 
(|4.68p . Thus Proposition is proved. □ 

4.3.3 The bulk contribution 

At this point we are left with the contribution coming from terms only containing Hao{B). Define: 

Xo{B,L) := Tr j fpoiz) { Pi{B)~go{Hoo{B) - z)-^goP2{B)~ga{H^{B) - z - uj)-^ga 

+ z-> z-Lj^dz. (4.102) 

We will compute voi(Ai.) ^^"^o(^' ^) show that it converges to {c^bctooKO). 
Define: 

Xo{B,L) :=Tr^ /^.^^(z) { Pi (5)50^5 (2)50/^2 (S)go5B(z + a;)go + z ^ z - w } dz. (4.103) 

Now we can prove the last technical result: 
Proposition 4.11. The following two double limits exist: 

'^i '-^ l™n 4{^o(i3, L) - Xo{B, L)}, (4.104) 

L^oo VoI(Al) b^o B 

■■= 1™ T^T(^-^ j.™n 4{^o(P, L) - Xo{0, L)}. (4.105) 
L^oo VoI(Al) b^o B 
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Moreover, the mapping sb defined in Theorem ] 1.1\ is differentiable at B = and 
lim {dsaLW) lim —L^^{dsXomL) 

= (71+^2 = {asdooKO) - - / {9BSs}B=o(x)dx. (4.106) 

Jn 

Proof. Let us start with (|4.104p . Using ((3?24)) one can show the foUowing identity: 
lim ^{Xo{B, L) - XoiB, L)} = ~Tr / fpoiz) 

■ { Pi(0)5o(i?oo(0) - z)-^ro(z)go/'2(0)5o(i?oo(0) - z - u;)-'go 
+ z— i-z — — Tr J fpoiz) 

■ { Pi(0).9o(i?oo(0) - z)-^goP2{0)go{H^{0) -z- uj)-^To{z + Lo)go 

+ z^z-uj^dz. (4.107) 

Then by integrating many times by parts, the integrand will become trace class, and we can get rid 
of the cut-ofF functions 50 and go since their removal will only contribute with a surface correction. 
Hence we can write: 

i^oo Voik) ^^^^'^ 

■ { Pi(0)(i/oo(0) - z)-^To{z)P2mH^iO) ~ z - w)-i 
+ z^z-a;}dz-^lkn^^^-^TrxA. IjMz) 

■ { Pi(0)(i/oo(0) - z)-iP2(0)(i?oo(0) - z ~ Lu)-'To{z + io) 
jrfz. (4.108) 



z — > z — a; 



Now one can prove (as we did for Foo) that the two operators defined above by integrals over Tui 
have jointly continuous integral kernels, whose diagonal values are Z^-periodic. It means that the 
limit exists and equals the integral of the kernels' diagonal value over the unit cube in M.^. 

Now let us continue with the proof of (|4.105[) . First, we can get rid of go because Pj{B) is 
local. Let us integrate by parts N times with respect to z, TV large. Then a typical term in the 
integrand defining Xo{B, L) will be: 

P^{B)S^^+'-''\z)goP2{B)S^B^'\z + oJ)go, k e {0,...,N}, 

where as before S^j^\z) has the integral kernel e^^'^°''^'^ ■'GAr(x, x'; z). 
This operator will have an integral kernel given by: 



Pi.x(S)e^^^''('''^)GAr+i_fc(x, y; z)go(y) 

•P2,y(S)e^^^''('''"''Gfc+i(y,x';z + ^)5o(x')dy. (4.109) 
We commute the momenta with the magnetic phases and obtain: 

f g.B^o(x,y)|p^^^(o) ^ _ y)}GAr+i_fe(x,y; z)go(y) 

g.B^o(y,x')|p^^^(0) + BA^iy - x')}Gfc+i(y,x'; z + Lo)go{^)dy. (4.110) 
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This integral is absolutely convergent and defines a continuous function in x and x' (we can see this 
from the regularity and exponential localization of GAr(x, x'; z) and its first order derivatives). In 
order to perform the trace of this operator we put x = x'. The two magnetic phases will disappear, 
thus we get: 

/ {Fi,x(0) + BAi(x - y)}Gjv+i-fc(x,y; z)go{y) 

{F2,y(0) + BA^iy - x)}Gfe+i(y,x; z + w).9o(x)dy. (4.111) 
The contribution to lims^o ^{^o{B, L) — Xo{0, L)} coming from this term will be: 

-Rl(x) (4.112) 
:=5o(x) / Ai(x-y)GAr+i_fc(x,y;z)5o(y)^2(y-x)Gfc+i(y,x;2+w)dy. 

Now we have to investigate the existence of the limit: 

1™ u^ / RL{x)dx. (4.113) 



L^oo VoI(Al) Jf, 

Let us first note that due to the exponential localization of G^'s (see (I3.1|) ) we have the following 
uniform estimate: 

sup / |Ai(x-y)| |GAr+i_fc(x,y;z)| |A2(y-x)| |Gfe+i(y,x;z+w)|dy 

< const •r*^ (4.114) 

If we look back at the definition of go, we see that it equals 1 on the complementary in of 
a boundary neighborhood like ^^(io) with to > 1 (see (l4.ip V Denote by xl the characteristic 
function of \ SL(2io)- Thus we have 

XL go = XL, dist{supp(l - go), supp(xL)} > toL". (4.115) 

Because of the uniform estimate (|4.114|) , the limit in (|4.113|) exists if and only if the following one 
exists: 

1™ TTTTTT / XL(x)i?L(x)dx, (4.116) 



L^oo VoI(Al) jaz. 

because the difference between integrands only gives a surface contribution. Let us now define: 
i?oo(x) (4.117) 
:= / Ai(x-y)GAr+i-fe(x,y;z)A2(y-x)Gfc+i(y,x;z + w)(iy. 

The difference between xlRl and xlRoo comes from the integration over the support of 1 — go- 
But due to (|4.115p and the exponential decay of Gfc's, this difference is of order e~^^'^/^'^\ thus 
will not contribute to the limit. Moreover, i?oo is Z'^-periodic, therefore we can write: 



lim sup (r) 



-M 



UA ^ / RL{^)dx~ / i?oo(x)dx 
VoI(Al) Jn 



0, (4.118) 



where M is some large enough positive number. Then the exponential decay of fpo will insure 
the convergence of the F^^-integrals, thus (|4.105p is proved. 

The last ingredient in the proof of Proposition 14.111 is the computation of c?b(Too(0) and the 
comparison with ai + (T2 • But the steps are very similar to those we have already done in order 
to compute cti and CT2. First, one integrates by parts many times with respect to z in order to 
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obtain a "nice" form for Fao- Second, using the magnetic perturbation theory one writes down a 
Taylor expansion in i? of (x) at i? = which only contains "regularized" terms and where we 
can interchange the expansion in B with the thermodynamic limit L oo. This strategy has been 
already used in [T3j for the Faraday effect (including the spin contribution, neglected here), and 
in [7] for generalized susceptibilities. □ 



5 Appendix: Uniform exponential decay 

The following proposition contains two key estimates which we are going to use throughout this 
paper. 

Proposition 5.1. Assume that z ^ C and dist{z, [0, oo)} — rj > 0. Then for any a G {1, 2, 3} we 

have 

sup \\[D^ + Baa] { (-«V + Ba)^ - z} ~' 1 1 

L>1 



< v/l^ + max{3fi(z),0}/r/2. (5.1) 
Moreover, there exists a constant C such that: 

sup sup {\z\)-^\\[Da + Baa]iHL{B)~z)-^\\<C/Tj. (5.2) 

L>1 sR(z)>0 

Proof. The estimate ()5.ip is an easy consequence of the following trivial identity, valid for every 

3 

Y,\\iDa+ Ba^] { (-zV + Ba)l - z}"' VI P (5.3) 
- (({(-zV + Bai)l - + 5R(^)||{(-*V + Bsi)l - z}-'^\\\ 

The estimate (|5.2p is a bit more involved. From (|5.ip we have that for every A > 1: 



C 

snpWiDa + Boc,] {{-iW + Bsi)l - tXy^ W < 
L>1 V A 

Since V is bounded we have 



sup ||F[(-iV + Ba)l - i\]'^\ \ < ^. 

L>1 A 

Choosing a Aq large enough and using the Neumann series in V for the resolvent we have 
sup {\\V{HLiB) - iXo)-'\\ + ||[(-zV + Ba)l - i\o\{Hl{B) - < 1/2. 

L>1 

Using the resolvent identity we obtain: 

sup ||[(-zV + Ba)l - iXo\{HL{B) - zY^W < C\z\/r]. (5.4) 

L>1 

Hence writing 

[Da + Baa][HL{B) - = [D^ + Ba„][(-zV + Bsi)l - iXo^^ 

■ [{-iV + Ba)!, - i\o]{Hl{B) - z)"i 

we obtain the result. □ 
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We will need a certain type of uniform exponential localization, stated in the proposition below. 
If Xq is some point in and a € M, then let e"^ "''"^ denote the multiplication operator with 
the function e^Vl'^^^^ol^+i. ]\fQ^g that multiplication with the exponential weight is a bounded 
operator if L < oo and leaves invariant the domain of Hl{B). 



Proposition 5.2. Fix xq E Al and dist{z, [0,oo)} = 77 > 0. Denote by r :— 5R(z). 
exists a 60 > and a constant C such that for every < S < Sq we have 



Then the 



sup sup sup 

L>1 rev. xoSA 



<c, 



(5.5) 



sup sup sup 

L>1 reRxoGA 



{(r)-i II [D,, + Ba^]e'^-^'>^^{HL{B) - z)-ie<-"°>T^ || } < C, (5.6) 



and 



{(r)-i ||e<--"°>T^f[(-iV + Ba)l + 1]{Hl{B) - z)-h^-''"^^j'^ < C. (5.7) 



sup sup sup 

L>1 rGRxoGA 

Proof. An heuristic explanation of (|5.5p is the following: if we apply the resolvent on a func- 
tion which is exponentially localized near Xq and decays like e~~^^^~^°\ then we do not loose 
exponential decay. Moreover, the bounds are uniform in z, L and the location of Xq. 
that For s G K, the well-known Combes-Thomas rotation [S] gives: 

3 

ga(-xo>(^^(5) - z)e-^<-''"> = Hl{B) - z + s^WjiDj + Baj] + sVi + sVa, 

where Wj, Vi, V2 are bounded functions, uniformly in L and Xq. 

Now put s = S/{r), and use (|5.2p . If 6 is small enough, we get that uniformly in L, xq and r 
we have 

3 

\\{sY,w,[D, + Ba,] + sV, + s'V^HHUB) - z)-'\\ < 1/2, 



which gives 



{HL{B)-z)-'e 
3 



{Hl{B)-z)-' 



s J2 {Dj + Baj) + sVi + s^V^, 



(5.8) 



This implies (|5.5p . and together with (|5.2p we also get 

Let us now concentrate ourselves on the last estimate (j5.7p . Up to a commutation, (|5.6p gives 

supsup sup |(r)-i • e'''^'^°^^[Da + Baa]{HL{B) - z)-^e^^'^'^°^^ I < C"- (5-9) 

L>1 reRxoGA J 

Thus again up to a commutation, (|5.7p follows if we can prove 

supsup sup |(r)-i [{-iV + Bs)']j + l]e^-'''>^^{HL{B)- zy^e'^--''"^^ \ <C. (5.10) 

L>1 rSRxoGA J 

But this estimate follows from (|5.8p and (|5.4p . 

□ 
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Corollary 5.3. Let A > Aq > 0. Then there exists c > such that 



sup sup sup 

L>1 X>Xo xoSA 



< 



const(Ao) 
A 



(5.11) 



Proof. We use the key estimate (|5.ip for the case when rj > X and 3?(z) = —A < 0. This gives us 
I + Ba„)[(-iV + Ba.)l + X]-^\ \ < const/\/A 

hence 

\\{D^+Baa,)[HLiB) + A]-i|| < const/ \/A. 

Now we proceed as in (|5.8p and we get the result. Finally, note that by repeating the argument 
of Proposition 15 . 21 we can obtain a uniform estimate in A, L and Xg: 



[{-tV + Ba)j, + X]{Hl{B) + X) 



< const. 



(5.12) 
□ 



Proposition 5.4. The operator [(— iV + Ba.)jj + X] ^ has an integral kernel K]^{x,x') which is 
jointly continuous away from the diagonal x = x', and obeys the estimate 



g- VA|x-x'| 

|i^,.(x,x')| < (5.13) 



for every x 7^ x' in . 



Proof. The argument is based on several well known results. First, one uses the Feynman-Kac- 
Ito representation for the kernel of the semi-group e~*(~*^+^'*)o, t > (see [S]) and obtains a 
diamagnetic inequality in A^: 

e-*(-™'^)-(x,x') <e*^°(x,x') < (4^t)-'/'e-^^, x,x'eAi. 
Second, we perform a Laplace transform and obtain the result. □ 



Proposition 5.5. All the results in this section are also valid if the operators are defined on the 
whole space M'^ (formally L ~ 00). 

Proof. The argument relies on various standard limiting and cut-off arguments, which are neces- 
sary because the exponential growing factors do not invariate the operator domain of Hao{B). The 
most important ingredient (uniformity in L > 1 of all our previous estimates) has been already 
proved. □ 
Acknowledgments. H.C. acknowledges support from the Danish F.N.U. grant Mathematical 
Physics. 
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